Note on proton-antiproton suppression in 200 AGeV Au-Au collisions 



Thorsten Renlfl and Kari J. EskolsQ 
Department of Physics, P.O. Box 35 FI-40014 University of Jyvaskyla, Finland and 
Helsinki Institute of Physics, P. O. Box 64 FI-00014, University of Helsinki, Finland 

We discuss the measured nuclear suppression of p + p production in 200 AGeV Au-Au collisions at 
RHIC within radiative energy loss. For the AKK set of fragmentation functions, proton production 
is dominated by giuons, giving rise to the expectation that the nuclear suppression for p + p should 
be stronger than for pions due to the stronger coupling of giuons to the quenching medium. Using 
a hydrodynamical description for the soft matter evolution, we show that this is indeed seen in the 
calculation. However, the expected suppression factors for pions and protons are sufficiently similar 
that a discrimination with present data is not possible. In the high px region above 6 GeV where the 
contributions of hydrodynamics and recombination to hadron production are negligible, the model 
calculation is in good agreement with the data on p + p suppression. 

PACS numbers: 25.75.-q,25.75.Gz 



The STAR collaboration has measured the nuclear sup- 
pression ratio Rcp (yield for — 12% central collisions 
divied by yield for 60 — 80% peripheral collisions) for 
7T+ + 7T~ and p + p in 200 AGeV Au-Au collisions [H, 
finding that the suppression pattern of both pions and 
protons is very similar in the 6+ GeV momentum region. 
This is to some degree surprising, as the AKK fragmenta- 
tion functions Q (which roughly describe the STAR data 
on proton production in p-p collisions indicate that 
~ 80 % of p+p production comes from gluon jets whereas 
only ~ 40 % of pion production is gluon-driven. This is 
not so in the older KKP set of fragmentation functions 
[H which consequently underpredict proton production 
at STAR. However, giuons as colour octet states are ex- 
pected to undergo stronger energy loss in medium than 
quarks, thus the dominance of gluon jets in hard pro- 
ton production should map into a stronger suppression 
of p + p production as compared to n + + n~ production. 

It becomes thus a quantitative question in the energy 
loss model to check how large this difference in suppres- 
sion should be and if present data are able to resolve it. 
A framework to study energy loss within full expansion 
dynamics of soft matter has been presented in [a] and 
extended to a 3-d hydrodynamical model @ in (7J. We 
utilize this framework to address the above question. We 
stress that given the scenario describing pionic Raa, the 
change to study the suppression of protons is only the 
trivial change of the fragmentation function, no further 
parameters are adjusted. Furthermore, since the AKK 
fragmentation function overpredicts p + p production in 
the high pt region in p-p collisions (and the fragmenta- 
tion function of q(q) — > p(p) is well constrained), thus 
overemphasizing the role of giuons, we stress that our 
results are an upper limit for the expected difference be- 
tween proton and pion suppression. 

Let us now discuss the treatment of partons propagating 



through the medium: our calculation follows the BDMPS 
formalism for radiative energy loss [1[ using quenching 
weights as introduced by Salgado and Wiedemann [9|, 
Hoj |. The probability density P(xq, yo) for finding a hard 
vertex at the transverse position ro = (xo, yo) and impact 
parameter b is given by the product of the nuclear profile 
functions as 



P (#0,2/o) 



r A (ro + b/2)r A (r - b/2) 
T AA (h) 
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where the thickness function is given in terms of 
Woods-Saxon the nuclear density p^(r, z) as Ta(i") = 
J dzpAiy, z). If we call the angle between outgoing par- 
ton and the reaction plane <fi, the path of a given parton 
through the medium £(r) is specified by (i"o,</>) and we 
can compute the energy loss probability P(AE) pat h for 
this path. We do this by evaluating the line integrals 

/>oo />oo 

Wc (r »= / df£q(0 and (qL)(r o ,0) = / dfq(£) 
Jo Jo 

along the path where we assume the relation 

q(£)=K-2- e 3/4 (£) (cosh p - sinh p cos a) (3) 

between the local transport coefficient q(£) (specifying 
the quenching power of the medium) , the energy density 
e and the local flow rapidity p with angle a between flow 
and parton trajectory [ill [l2j]. Here oj c is the charac- 
teristic gluon frequency, setting the scale of the energy 
loss probability distribution, and (qL) is a measure of the 
path- length weighted by the local quenching power. We 
view the parameter K as a tool to account for the uncer- 
tainty in the selection of a s and possible non-perturbative 
effects increasing the quenching power of the medium (see 
discussion in Q) and adjust it such that pionic Raa for 
central Au-Au collisions is described. 
Using the numerical results of ||, we obtain 
P(AE;oj c , R) P ath for uj c and R = 2u1/{qL) as a func- 
tion of jet production vertex and the angle <fi. 
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FIG. 1: Raa for pions as measured by PHENIX 13|, scaled 
Rcp for pions measured by STAR li] (in good agreement 
with the PHENIX Raa), scaled Rcp for protons measured by 
STAR jl( and calculated Raa for pions and protons, dashed 
lines indicate the regime where recombination is expected to 
become an important effect [13 ]. 



FIG. 2: Calculated Raa for quarks and gluons and pro- 
tons/antiprotons originating from quarks and gluons using 
the AKK fragmentation functions as well as the final result 
averaged over the relative contributions of quarks and gluons 
in the pQCD calculation. 



From the energy loss distribution given a single path, we 
can define the averaged energy loss probability distribu- 
tion P(AE)) Taa as 



(P(AE)) Ta , 
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dx Q / dy P{x ,y )P(AE) 

pa 



th- 



(4) 

The energy loss probability P(AE) pat h is derived in the 
limit of infinite parton energy [9] . In order to account for 
the finite energy of the partons we truncate (P(AE))t aa 
at AE = E- ]ct and add S(A—Ej ei ) JJ° deP(e) to the trun- 
cated distribution to ensure proper normalization. The 
physics meaning of this correction is that we consider all 
partons as absorbed whose energy loss is formally larger 
than their initial energy. We calculate the momentum 
spectrum of hard partons in leading order perturbative 
QCD (LO pQCD) (explicit expressions are given in [5| 
and references therein). The medium- modified perturba- 
tive production of hadrons can then be computed from 
the expression 
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with Df°^ h (z, fj,p) the fragmentation function with mo- 
mentum fraction z at scale /j, 2 f Q , and from this we com- 
pute the nuclear modification factor Raa as 



Raa{pt,v) 



dN'XJdprdy 
T AA (h)daPP/dp T dy' 



(6) 



In Fig. [T] we compare with the STAR data. Since our 
model, making use of thermodynamical quantities to de- 
scribe energy loss induced by the soft medium is not 
suited to describe 60 — 80% peripheral collisions where 
a dilute, predominantly hadronic, medium is expected 



to exist, we calculate Raa and scale the experimental 
data. We show that the scaled pionic Rcp as obtained 
by STAR 1] is in perfect agreement with the pionic Raa 
as obtained by PHENIX pj|]. This gives us confidence 
that the comparison of protonic Raa with scaled Rcp is 
not grossly unreasonable. As apparent from the results, 
the similarity of the pionic and protonic suppression is 
expected, the difference between pion and proton Raa is 
less than the error on the pion Raa over a large kinematic 
range. 

Let us briefly discuss the origin of this outcome: While 
for an infinite energy parton the average energy loss scales 
linearly with the coupling strength to the medium (and 
hence with the color charge), Raa for finite pt cannot 
exhibit a linear scaling downward for the simple reason 
that it is bounded from below by zero. In other words, 
once the medium quenching is strong enough such that 
essentially all partons from a spacetime region are ab- 
sorbed, increasing the quenching power of this region 
further cannot alter Raa any more and consequently 
observed hadrons in such a situation are emitted from 
the medium surface (cf. the discussion in [TEl]). In the 
present calculation, we likewise observe this saturation 
effect (albeit somewhat weaker, cf. Q), predominantly 
for gluons. 

In Fig. [2] we show Raa for quarks and gluons as well 
as Raa for protons produced exclusively from quarks 
and gluons. Since the fraction of protons produced from 
quarks is small when using the AKK fragmentation func- 
tions, the final result for protons is very close to the gluon 
Raa at all accessible px (whereas the pion Raa is be- 
tween the two partonic curves, showing the transition 
from a gluon-dominated parton spectrum at low pt to 
a quark-dominated spectrum at higher pt)- The gluonic 
Raa is clearly in a saturated regime and we checked that 
the spread between quark and gluon Raa is reduced by a 
factor two as compared to a medium with less quenching 



3 



power in which Raa for pions is close to 0.5. It is also 
apparent that while the functional form of the fragmenta- 
tion function has some influence on the results at lower 
Pt, Raa can mainly be understood from the partonic 
suppression and the px dependent ratio of quarks and 
gluons contributing to the production of a given hadron. 
Thus, the similarity of Raa for protons and pions in this 
calculation has to be attributed to the fact that gluon 
emergy loss appears in the saturated regime. 
Once again, we stress that we consider the calculation to 
be an upper limit for the expected difference. For the 
most part, the STAR data on proton suppression are in 
a regime where hydrodynamics or recombination [14| is 



the dominant contribution to hadron production, in the 
momentum region above 6 GeV where fragmentation is 
expected to dominate hadronization, the calculation is in 
agreement with the data. 
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